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The study of jet production has become one of the precision tests of QCD. Next-to-leading order 
(NLO) computations have been successfully confronted with experimental data. For example, jet pro¬ 
duction at the TeVatron ||I|] provides a very stringent test of available QCD Monte Carlo codes and gives 
very valuable information on the behaviour of the gluon distribution in the proton at large x. 

On the other hand, jets have never been measured in a nuclear environment. There are prospects 
to measure high-p^ particles at RHIC, but only at the LHC the opportunity of detecting and studying 
jets with Et ~ 100 GeV will appear. Concretely, such measurements in pA collisions would allow 
to settle some uncertainties like the validity of collinear factorization [§], the modifications of nucleon 
parton density functions (pdf) inside nuclei |^, or the existence and amount of energy loss of fast partons 
inside cold nuclear matter |0]. Such issues have to become well determined before any claim, based on 
jet studies, on the existence of new physics in nucleus-nucleus collisions (hopefully QGP formation), 
can be considered conclusive. Furthermore, the study of high-p^ particles and jets can help not only to 
disentangle the existence of such new state of matter but also to characterize its properties, see e.g. [^. 

For our computations we will use the Monte Carlo code at NLO of adapted to include 

isospin effects and modifications of nucleon pdf inside nuclei. This code is based on the subtraction 
method to cancel the infrared singularities between real and virtual contributions; full explanations and 
a list of available codes, as well as detailed discussions on theoretical uncertainties in nucleon-nucleon 
collisions, can be found in [^]. The accuracy of our computations is limited by CPU time, but it can be 
estimated to be: 

• For the transverse energy distributions, 2 % for the lowest and 15 % for the highest i?r-bins. 

• For the pseudorapidity distributions, 3 %. 

• For the dijet distributions of the angle between the two hardest jets, 20 % for the least populated 
and 3 % for the most populated bins. 

The results will be presented in the LHC lab frame, i.e. that in which the experiments work. For example, 
for pPb collisions this means a 7 TeV proton beam against a 2.75 TeV Pb beam, see the Section on the 
experimental parameters in pA collisions. All the energies will be given per nucleon and, in order to 
compare with the pp case, all cross sections will be presented per nucleon-nucleon pair, i.e. divided by 
AB. 

Unless explicitly stated, we will use as nucleon pdf the MRST98 central gluon distribution 0 
modified inside nuclei using fhe EKS98 parameferizafions [jl^, a facforizafion scale equal fo fhe renor- 
malizafion scale /r = //j? = fiR = Et /2 (wifh Et the total transverse energy of all the jets in the gen¬ 
erated event), and for jet reconstruction we will employ the fcr-clustering algorithm [ |n| ] with D = 1, 
which is more sound on theoretical grounds than the cone algorithm [jT2|]. The kinematical regions we 
are going to consider are the following: 

• {rjil < 2.5, with rji the pseudorapidity of the jet; this corresponds to the acceptance of the central 
part of the CMS detector. 

• Eti > 20 GeV in the pseudorapidity distributions, with Eti the transverse energy of the jet; this 
will ensure the validity of perturbative QCD. 

• Eti > 20 GeV and Et 2 > 15 GeV for the (/)-distributions, with Eti {Eti) the transverse energy 
of the hardest (next-to-hardest) jet entering the CMS acceptance, and cj) the angle between these 
two jets. 


Please note that, even in the absence of nuclear effects, the 77 *-distributions may be asymmetric with 
respect to r/j = 0 due to the fact that we are giving our results in the lab frame, not in the center-of-mass 
frame. 


Concerning the dijet distributions in the angle between the two hardest jets, asymmetric cuts have 
been imposed to avoid the generation of large logarithms in certain points of the phase space, see [^. 
Also, the results near 0 = tt are not reliable 0 , as they require an all-order resummation not available 
in the code we are using. 


Some call for caution, which will be even more important for AB collisions (see the Section on Jet 
and Dijet Rates in AB Collisions [ 131), about the possible comparison of our results with the experimental 
situation has to be stressed at this point. There exist several physical effects that are not included in 
our computations. For example, the so-called underlying event (i.e. the soft particle production that 
coexists with the hard process) is not considered here; it may cause difficulties with jet reconstruction 
and increase the uncertainties in jet-definition algorithms. This effect has been considered in antiproton- 
proton collisions at the TeVatron [ 14 ] but its estimation, which relies on our limited knowledge of soft 
multiparticle production, is model dependent. In collisions involving nuclei the situation is even worse, as 
our knowledge is more limited. The only way to take this into account would be to use available Monte 
Carlo simulators for the full pA event [15], suitably generalized to introduce NLO QCD calculations. 
Another process not included in our computations is multiple hard parton collisions [p^. They may be 
classified info fwo classes: 


a) Disconnecfed collisions, i.e. processes in which in fhe same even! fhere are more fhan one in¬ 
dependent parton-parton collisions producing a pair of high-px jets each. This aspect has also 
been studied at the TeVatron [ JT/I ] for hard collisions coming from one nucleon-nucleon collision, 
but its quantitative explanation and the extension to collisions involving nuclei is model depen¬ 
dent. Simple estimates of the influence of disconnecfed collisions on jef production in pA col¬ 
lisions may be obfained by computing fhe number (n) of nucleon-nucleon collisions involved in 
fhe production of jefs wifh Eti greater than a given Eto- In the Glauber model [ |I^ one obtains: 
(n)(5, Eto) = ATA{b)cF{ETo)/CFpA{b, Etq), with b the impact parameter, T^( 6 ) the nuclear pro¬ 
file function normalized fo unify, a {Etq) the cross section for production of jets with Eti greater 
than Etq in pp collisions, and apA{b,ETo) = 1 — [1 — T^( 6 )(t(£^to)]^- Taking a{ETo) = 120, 
0.3 and 0.02 pb as representative values in pPb collisions at 8.8 TeV for Etq = 20, 100 and 200 
GeV respectively (see results in Fig. ^ below), the number of nucleon-nucleon collisions involved 
turns out to be 1.0 for all Etq, both for minimum bias collisions (i.e. integrating numerator and 
denominator in apA{b, Eto) between 6 = 0 and 00 ) and for central collisions (e.g. integrating 
between 6 = 0 and 1 fm). So, in pA collisions at LHC energies the contribution of multiple 
hard scattering coming from different nucleon-nucleon collisions seems to be negligible for the 
transverse energies of the jets considered in this study. 

b) Rescatterings, i.e. processes in which a given high-p-r parton may undergo several hard collisions 
before hadronizing into a jet, mimicking a single process of higher order in the perturbative QCD 
expansion. Quantitative descriptions of the resulting modification of the jet transverse energy 
spectrum are model dependent. However, the effect at LHC is in general expected to be very 
small, and almost negligible at the transverse energies considered in this paper (see the Section on 
Cronin effect in proton-nucleus collisions: a survey of theoretical models [|I^]). 

As a last comment, no centrality dependence has been studied. It is not clear how to implement such 
a dependence in theoretical computations, because not only the number of binary nucleon-nucleon col¬ 
lisions should change with centrality, but also the modification of nucleon pdf inside nuclei may be 
centrality-dependent. 





1. do jdEo and do / dr] for jets at large Et 

1.1 Uncertainties 

Four sources of theoretical uncertainties in pp and pPb collisions are examined in Figs. 0 and ^ First, it 
can be seen that varying the choice of the scale between Et/4: and Et gives differences of order ±10 %, 
the smaller scale giving the larger results. Second, the variation due to a different choice of nucleon 
pdf has been examined by using CTEQ5M as nucleon pdf [^], which produces larger results than the 
default choice MRST98 central gluon 0, the variation being of order 5 %. Third, the change due to 
isospin effects (obtained from the comparison of pp and pPb without any modifications of nucleon pdf 
inside nuclei at the same energy per nucleon) is very small as one would expect from the dominance of 
the gluon-gluon production channel; the effect of modifications of nucleon pdf inside nuclei, estimated 
by using EKS98 [ |T^ ] nuclear corrections, is also small but produces an extra asymmetry (an excess in 
the region r]i < 0) in the pseudorapidity distributions for pA collisions at reduced energy (e.g. pPb at 
5.5 TeV), which disappears at maximum energy (e.g. pPb at 8.8 TeV), see Eig. ^ (bottom-left) and Eig. 
^ (bottom-right), due to the the asymmetry of momentum of projectile and target in the EHC lab frame. 
Einally, the choice of jet-finding algorithm produces some differences: the cone algorithm [ |T^ with 
R = 0.7 gives results slightly smaller (~ 1 2 %) than the /sT-dustering algorithm [|^] with D = 1, 

while cone with R = 1 gives results ~ 15 % larger than our default choice. 

As a last comment, let us discuss the ratio of cross sections evaluated at NEO over those computed 
at leading order (EO), the so-called AT-factor. In Eig. || this ratio has been examined for different 
collisions and energies, and varying the different sources of uncertainties. Eor our default option the 
ratio results quite constant with the transverse energy or pseudorapidity of the jet and turns out to be 
around 1.2 for the energies examined. The dependence on the renormalization/factorization scale results 
in a ratio ~ 1.35 for ^ = Et and ~ 1.15 for /r = Et/4:. The cone algorithm with R = 0.7 gives 
results very similar to those obtained with the /cr-clustering algorithm with D = 1, while cone with 
R = I produces a ratio ~ 1.45 (note that at EO the choice of jet-finding algorithm has no influence on 
the results). Einally, neither a variation of nucleon pdf, nor isospin effects and modifications of nucleon 
pdf inside nuclei, have a sizeable effect. 

1.2 Results 

In Table 1 the number of expected events with at least one jet with a given Eti > 20 GeV and \rii\ < 2.5 
(or with two jets (1, 2) with Eti > 20 GeV, Et 2 > 15 GeV and |7/i 2 1 < 2.5 for the dijet (/)-distributions), 
per ]ih and pair of colliding nucleons is shown for different collisions and possible luminosities. Erom 
this Table and using the Eigures it is possible to know the number of expected events with a given 
kinematical variable. Eor example, examining the solid line in Eig. ^ (upper-left) one can expect, within 
the pseudorapidity region we have considered, the following numbers of jets in pp collisions at 14 TeV 
for a luminosity of 10^^ cm“^s“^: 10^° jets with Eti ~ 70 GeV (corresponding to a cross section of 1 
^b/(AB)), and 10® jets with Eti ~ 380 GeV (corresponding to a cross section of 10 ^ /rb/(AB)). 

Booking at the results given in Eig. it becomes evident that the study of samples of ~ 10^ jets 
should be feasible, from a theoretical point of view, up to a transverse energy ETi ~ 325 GeV with a run 
of 1 month at the considered luminosity: indeed, looking at Table 1, 10^ jets for pPb at 8.8 TeV would 
correspond to a cross section of 4.8 ■ 10“® /rb/(AB), which in Eig. ^(upper-right) cuts the solid curve at 
ETi ~ 325 GeV. 

Some conclusions can be obtained at this point: Eirst, as the influence of modifications of nucleon 
pdf inside nuclei is not large, an extensive systematic study of the A-dependence of these cross sections 
does not seem to be the most urgent need. Second, noting the asymmetry in the pseudorapidity distribu¬ 
tions in pPb collisions at the reduced energy of 5.5 TeV (and to a lesser extent in pAr at 6.3 TeV) when 
EKS98 corrections are implemented, a pA run at reduced proton energies might offer a possibility to test 
the influence of the modification of nucleon pdf inside nuclei; a study of A-rj-distributions in different 







Fig. 1: Plots on the left: Scale dependence of jet cross sections (/r = EtI2: solid lines; /r = Erji: dashed lines; jj, = Et'- 
dotted lines) versus transverse energy of the jet (for \rji\ < 2.5, upper plot) and pseudorapidity of the jet (for En > 20 GeV, 
lower plot). Plots on the right: Nucleon pdf dependence of jet cross sections (MRST98 central gluon: solid lines; CTEQ5M: 
dashed lines) versus transverse energy of the jet (for \r]i\ < 2.5, upper plot) and pseudorapidity of the jet (for En > 20 GeV, 
lower plot). In each plot, results for pPb at 8.8 TeV (upper lines) and pp at 14 TeV (lower lines, multiplied by 0.1 in the En-piot 
and by 0.25 in the r^i-plot) are shown. Unless otherwise stated default options are used, see text. 


r)j-slices going from the backward (A) to the forward (p) region would scan different regions in x of the 
ratio of the gluon distribution inside the nucleus over that in a free nucleon, which is badly under control 
with the presently available experimental information. Nevertheless, this would require a detailed com¬ 
parison with pp results at the same energy in order to disentangle possible biases (e.g. detector effects); 
in case these pp results were not available experimentally, the uncertainty due to the scale dependence 
(which in our computations at reduced energies is as large as the asymmetry) in the extrapolation from 
14 TeV to the energy of the pA collision should be reduced as much as possible for this effect to be 
useful. 

2. High-E’r dijet momentum imbalance 

Dijet distributions of the angle between the two hardest jets offer a possibility to test the perturbative 
expansion in QCD. At lowest order (LO) in collinear factorization the two jets are produced back- 
to-back, so any deviation from a simple peak at = vr in our results is a signal of NLO corrections. 
But please remember that the results near 4> = tt are not reliable [0], as there NLO corrections become 
negative and larger than LO, so this region requires an all-order resummation. 







Fig. 2: Plots on the left. Isospin and nuclear pdf dependence of jet cross sections (pp results: solid lines; pA results without 
modification of nucleon pdf inside nuclei: dashed lines; pA results with EKS98 modification of nucleon pdf inside nuclei: 
dotted lines) versus transverse energy of the jet (for |77i| < 2.5, upper plot) and pseudorapidity of the jet (for En > 20 GeV, 
lower plot). In each plot, results for pp and pPh at 5.5 TeV (upper lines) and pp and pAr at 6.3 TeV (lower lines, multiplied 
by 0.1 in the i?Ti-plot and by 0.25 in the T^i-plot) are shown. Plots on the righf. Dependence of jet cross sections on the jet 
reconstruction algorithm (fer-algorithm with D = 1\ solid lines; cone algorithm with R = 0.7: dashed lines; cone algorithm 
with R = 1\ dotted lines) versus transverse energy of the jet (for \r]i\ < 2.5, upper plot) and pseudorapidity of the jet (for 
En > 20 GeV, lower plot). In each plot, results for pPb at 8.8 TeV (upper lines) and pp at 14 TeV (lower lines, multiplied by 
0.1 in the ifCTi-plot and by 0.25 in the j^i-plot) are shown. Unless otherwise stated default options are used, see text. 


2.1 Uncertainties 

In Fig. ^ the same uncertainties examined for transverse energy and pseudorapidity distributions are 
studied in the dijet cross sections versus angle between the two hardest jets. First, the smaller the scale, 
the largest the results, so the results for p- = Et become up to ~ 50 % smaller than for /r = Et/^- 
Second, the choice of nucleon pdf seems to have a small effect (less than 10 %), as we have tested by 
using CTEQ5M instead of MRST98 central gluon. Third, the variation due to isospin effects (obtained 
from the comparison of pp and pPb without any modifications of nucleon pdf inside nuclei at the same 
energy per nucleon) is again negligible, and the effect of modifications of nucleon pdf inside nuclei, 
estimated by using EKS98 [10| nuclear corrections, is also small (~ 3 %). Einally, the choice of jet¬ 
finding algorifhm produces some differences which can be of order 20 % for our fhree choices of cone 
wifh R= 1, cone wifh R = 0.7 and fhe defaulf /cT-dustering algorifhm wifh D = 1. 
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Fig. 3: Ratio of cross section computed at NLO over those computed at LO versus the transverse energy of the jet (for \rii \ < 2.5, 
plots on the left) and the pseudorapidity of the jet (for ETi > 20 GeV, plots on the right). From top to bottom: scale dependence 
(/i = EtI2: solid lines; pt = EtI^’. dashed lines; pt = Et- dotted lines); dependence on the jet reconstruction algorithm 
(fer-algorithm with D = 1: solid lines; cone algorithm with R = 0.7: dashed lines; cone algorithm with R = 1: dotted lines); 
nucleon pdf dependence (MRST98 central gluon: solid lines; CTEQ5M: dashed lines); and isospin and nuclear pdf dependence 
(pp results: solid lines; pA results without modification of nucleon pdf inside nuclei: dashed lines; pA results with EKS98 
modification of nucleon pdf inside nuclei: dotted lines). Six upper plots: results for pp collisions at 14, 6.3 (multiplied by 0.5) 
and 5.5 (multiplied by 0.2) TeV are shown. Plots at the bottom: results for pp and pAr collisions at 6.3 TeV (multiplied by 0.5), 
and for pp and pPb collisions at 5.5 TeV, are given. Unless otherwise stated default options are used, see text. 


2.2 Results 

Results for pp and pA collisions at different energies are presented in Fig. From this Figure and Table 
1, the dijet momentum imbalance should be clearly measurable, provided jet reconstruction is possi¬ 
ble in the nuclear environment and no other physics contribution, like the underlying event or multiple 
hard parton scattering, interferes, spoiling the comparison with experimental data of the theoretical pre¬ 
dictions we present. Again, only extensive studies using Monte Carlo simulators including full event 
reconstruction will be able to clarify whether such measurements are feasible or not, see the Section on 
Jet Detection at CMS. 































Table 1: Luminosities and expeeted number of events with at least one jet with a given Eti > 20 GeV 
and |? 7 j| < 2.5 (or with two jets (1, 2) with Eti > 20 GeV, Et 2 > 15 GeV and |r/i, 2 | < 2.5 for the dijet 
(^-distributions), per ;ub/(AB) in one month (10® s), for different eollisions. 


Collision 

Ecra per nueleon (TeV) 

C (em ^s ^) 

Number of events per month per /ib/(AB) 

PP 

14 

I(F 

To™ 

PP 

14 

3 • 10®® 

3-10® 

pAr 

9.4 

4 • 10®® 

1.6 • 10« 

pAr 

9.4 

1.2 • lO^y 

4.8 • 10® 

pPb 

8.8 


2.1 • 10'' 
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Fig. 4: Energy dependence of jet cross sections versus transverse energy of the jet (for \rii\ < 2.5, upper plots) and pseudora¬ 
pidity of the jet (for Eti > 20 GeV, lower plots). Plots on the left, pp collisions at 14 (solid lines), 6.3 (dashed lines) and 5.5 
TeV (dotted lines). Plots in the middle: pAr collisions at 9.4 (solid lines) and 6.3 TeV (dashed lines). Plots on the right: pPb 
collisions at 8.8 (solid lines) and 5.5 TeV (dashed lines). Default options are used, see text. 
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Fig. 5: Uncertainties in dijet cross sections versus angle between the two hardest jets for Eti > 20 GeV, Et 2 > 15 GeV and 
|7?i|, \r} 2 \ < 2.5. Upper left plot, scale dependence in pp collisions at 14 TeV with the same line convention as in Fig. I^eft. 
Upper right plot, dependence on the jet reconstruction algorithm in pp collisions at 14 TeV with the same line convention as in 
Fig. ^ight. Lower left plot: nucleon pdf dependence in pp collisions at 14 TeV with the same line convention as in Fig. |^ight. 
Lower right plot: isospin and nuclear pdf dependence in pp and pPb collisions at 5.5 TeV with the same line convention as in 
Fig. ^eft. Unless otherwise stated default options are used, see text. 
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Fig. 6: Energy dependence of dijet cross sections versus angle between the two hardest jets for Eti > 20 GeV, Et 2 > 15 
GeV and \t]i\, \ri 2 \ < 2.5. Upper plot: results in pp collisions with the same line convention as in Fig. ^eft. Medium plot: 
results in pAr collisions with the same line convention as in Fig. ^niddle. Lower plot: results in pPb collisions with the same 
line convention as in Fig. 0-ight. Default options are used, see text. 












